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It has been shown that tribimaximal mixing can be obtained by some particular breaking 
pattern of the symmetry, wherein the extra triplet Higgs scalars pick up certain fixed vacuum 
O expectation value (VEV) alignments. We have performed a detailed analysis of the different 
possible neutrino mass matrices within the framework of the model. We take into account 
all possible singlet and triplet Higgs scalars which leave the Lagrangian invariant under A^. We 
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^ break A4, spontaneously, allowing the Higgs to take any VEV in general. We show that the 



neutrino mixing matrix deviates from tribimaximal, both due to the presence of the extra Higgs 
singlets, as well as from the deviation of the triplet Higgs VEV from its desired alignment, taken 
previously. We solve the eigenvalue problem for a variety of these illustrative cases and identify 
the ones where one obtains exact tribimaximal mixing. All such cases require fine-tuning. We 
show which neutrino mass matrices would be strongly disfavored by the current neutrino data. 
Finally, we study in detail the phenomenology of the remaining viable mass matrices and establish 
the deviation of the neutrino mixing from tribimaximal, both analytically as well as numerically. 



1 Introduction 



Neutrinos have provided us with a window to physics beyond the Standard Model. A plethora 
of striking experimental results have propelled us to a juncture where we already know a great 
deal about the basic structure of the neutrino mixing matrix. Results from KamLAND [1] and 
solar neutrino experiments [2] can be best explained if sin^ 9i2 — 0.32, while atmospheric neutrino 
experiments [3] and results from K2K [4] and MINOS [5] pick sin^ 623 — 0.5 as the best-fit solution. 
So far there has been no evidence for ^13 driven oscillations and hence ^13 is currently consistent 
with zero, with a 3a upper bound of sin^ ^^13 < 0.05 [6, 7]. These results prod us to believe that 
the mixing matrix in the lepton sector could be of the tribimaximal (TBM) mixing form, first 
proposed by Harrison, Perkins, and Scott [8, 9], 
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It has been well known that for TBM mixing to exist, the neutrino mass matrix should be of the 
form 



/A B ^ \ 

B l(A + B + D) l(A + B-D) 
\B l(A + B-D) l(A + B + D)J 



(2) 



where A = |(2mi + 777.2), B = ^{1712 — mi) and D = rris, where 7771, 7772 and 7773 are the neutrino 
masses. The current neutrino data already give very good measurement of mass squared differ- 
ences, while the best limit on the absolute neutrino mass scale comes from cosmological data. We 
show in Table 1 the allowed neutrino oscillation parameters within their 3a ranges. 
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sin^ 6*12 


sin^ 023 


sin^ 013 


(7.1 - 8.3) X 10-5 eV^ 


(2.0 - 2.8) X 10-3 eV^ 


0.26 - 0.40 


0.34 - 0.67 


< 0.05 



Table 1: The 3a allowed intervals (1 dof) for the three-flavor neutrino oscillation parameters from 
global data including solar, atmospheric, reactor (KamLAND and CHOOZ) and accelerator (K2K 
and MINOS) experiments. 



Maximal 623 and zero ^13 can be easily obtained if Ai^, possesses /7 — r exchange symmetry [10] 
or the — Lr symmetry [11]. However, the solar mixing angle sin^ 6*12 is not so easily predicted to 
be exactly 1/3, as required for exact TBM mixing. While the current data points towards Eq. (1), 
better precision on the mixing angles are needed to really test the TBM ansatz and deviation from 
TBM [12] . The mixing angles 613 will be probed in the next generation long baseline and reactor 
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experiments [13], while deviation of ^23 from maximality could be done in atmospheric neutrino 
experiments [14]. However, it is extremely crucial that one makes a very accurate measurement 
of 9i2 in the future [15] in order to confirm TBM mixing. 

The challenge for model builders lies in explaining all features of the lepton mixing matrix 
together with the mass pattern of the neutrinos given in Table 1. If indeed the neutrinos have 
TBM mixing, then one should be able to naturally generate the mass matrix given in Eq. (2). 
Various symmetry groups explaining the flavor structure of the leptons have been invoked in the 
literature in order to accommodate the neutrino masses and mixing along with the charged leptons. 
In particular, the study of the non-Abelian discrete symmetry group A4 has received considerable 
interest in the recent past [16, 17, 18, 19, 20, 21, 22, 23]. This group has been shown to successfully 
reproduce the TBM form of the neutrino mixing matrix, in the basis where the charged lepton 
mass matrix is diagonal [19] . However, the authors of [19] work in a very special framework where 
only one of the three possible singlet Higgs under A4 is considered and for the triplet Higgs 
which contributes to the neutrino mass matrix, a particular vacuum alignment is taken^. In this 
framework, the mixing matrix emerges as apparently independent of the Yukawa couplings, the 
VEVs and the scale of A4 breaking. Only the mass eigenstates depend on them. 

In this paper, we consider the most general scenario with all possible Higgs scalars that can 
be accommodated within this A4 model. We give both approximate analytical solutions of the 
expected phenomenology as well as exact numerical results. We expound the conditions on VEVs 
and Yukawas needed for obtaining exact TBM mixing and show how the neutrino masses in those 
cases severely constrain them. We begin by showing that in the model considered in [19], one 
gets TBM mixing simply through the alignment of the A4 triplet Higgs vev. A concerted effort 
has been made in the literature to explain naturally this particular vacuum alignment needed for 
TBM mixing. The mass squared differences can be obtained if one has the vacuum expectation 
value (VEV) of an additional singlet Higgs. We show that to get the correct Amgi and Am|j^, the 
product of the VEV and Yukawa of this singlet is determined almost completely by the VEV and 
Yukawa of the triplet. This emerges as a further undesirable feature of the model, and one would 
need further explanation for this additional "relative ahgnment" between the product of VEVs 
and Yukawas of the triplet and the singlet. 

We allow for presence of an additional singlet Higgs, construct the neutrino mass matrices and 
study their phenomenology. We take one, two and three singlet Higgs at a time and check which 
combinations produce viable neutrino mass matrices. In particular, we check which ones would 
produce TBM mixing and under what conditions. We obtain a few combinations in addition to 
the one considered in [19] which produce TBM mixing. We perturb these mass matrices by slightly 
shifting the VEVs and/or the Yukawas and study the deviation from TBM and the corresponding 
effect on the neutrino masses. We note that if just one singlet Higgs is allowed, the model of 
[19] is the only viable model. We further show that in the simplest version of this model, one 
necessarily gets the normal mass hierarchy^. Inverted hierarchy can be possible if we have at least 
two or all three Higgs scalars with nonzero VEVs. We give predictions for the sum of the absolute 
neutrino masses (mi), the effective mass that will be observed in neutrinoless double beta decay 

^We give a brief overview of the A4 model in section 2. 

^We call this the neutrino mass hierarchy, though what we mean is the ordering of the neutrino mass states. 
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Table 2: List of fermion and scalar fields used in this model. Two lower rows list additional singlets 
considered in the present work. In section 4, we also allow for a different VEV ahgnment for 05. 



experiments ((mge)) and for tritium beta decay experiments. Finally, we perturb the VEV 
alignment for the triplet Higgs and study the deviation from TBM mixing. 

In section 2 we give a brief overview of the model considered. In section 3 we begin with 
detailed phenomenological analysis of the case where there is just one singlet Higgs under ^44. We 
next increase the number of contributing singlet Higgs , give analytical and numerical results. In 
section 4 we study the impact of the misalignment of the VEVs of the triplet Higgs. We end in 
section 5 with our conclusions. 



2 Overview of the Model 

Alternating group is a group of even permutations of n objects. It is a subgroup of the 
permutation group Sn and has y elements. The non-Abelian group A4, is the first alternating 
group which is not a direct product of cyclic groups, and is isomorphic to the tetrahedral group 
Td- The group A^ has 12 elements, which can be written in terms of the generators of the group 
S and T. The generators of satisfy the relation 

There are three one-dimensional irreducible representations of the group ^4 denoted as 

(4) 
(5) 
(6) 
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It is easy to check that there is no two-dimensional irreducible representation of this group. The 
three-dimensional unitary representations of T and S are 

/100\ ^ / -1 2 2 \ 

T = 1 5 = - 2 -1 2 . (7) 

\0 u J ^ \ 2 2 -1 J 

where T has been chosen to be diagonal. We refer the readers to [19] for a quick review of the A4 
group. Here we give multiplication rules for the singlet and triplet representations, for the sake of 
completeness. These multiplication rules correspond to the specific basis of two generators S, T of 
A4. We have 

1x1 = 1, l'xl" = l 3x3 = 3 + 3a + 1 + 1' + 1" . (8) 

For two triplets 

a = (01,02,03), (61, 62, 63) (9) 

one can write 

I = (ab) = (oifei + 02^3 + 03^2) (10) 
l' = (o6)' = (0353 + 0162 + 02^1) (11) 
l" = (o6)" = (0262 + 0163 + 0361). (12) 

Note that while 1 remains invariant under the exchange of the second and third elements of o and 
6, 1' is symmetric under the exchange of the first and second elements while 1" is symmetric under 
the exchange of the first and third elements. 

3 = (06)5 = ^ |^2oi6i - 0263 - 0362, 20363 - 0162 - 0261, 20262 - 0163 - 0361^ (13) 
3^ = (06)^ = ^("0263-0362,0162-0261,0163-0361). (14) 
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We will be concerned with only 3 and we can see that the first element here has 2-3 exchange 
symmetry,the second element has 1-2 exchange symmetry, while the third element has 1-3 exchange 
symmetry. ^ 

The full particle content of the model we consider is shown in Table 2. There are five SU (2) (8) 
C/y(l) Higgs singlets, three (^, ^' and of which are singlets under A4 and two {(px and (ps) 
of which transform as triplets. The standard model lepton doublets are assigned to the triplet 
representation of A4, while the right handed charged leptons e'^, /x'^ and r"^ are assumed to belong 
to the 1, 1" and 1' representation respectively. The standard Higgs doublets /i„ and hd remain 
invariant under A4. The form of the A4 invariant Yukawa part of the Lagrangian is^ 

jCy = Vee'iM + y,f^%(f>Tiy + Vrr'iM" + + W + W + Xbicf>sU) + h.c. + ... (15) 



^There are no 3a terms in the Lagrangian since it is antisymmetric and hence cannot be used for the neutrino 
mass matrix. 

"^Wc assume that ips docs not couple to charged leptons and (f>T docs not contribute to the Majorana mass 
matrix. These two additional features can be obtained from extra-dimensional realization of the model, or from 
extra abelian symmetries [19]. 
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where, following [19] we have used the compact notation, yf-e^^cpxl) = yeG'^{4>Ti)hd/ A, Xa^ill) = 
Xaiilhylhu) / and so on, and A is the cut-off scale of the theory. After the spontaneous breaking 
of Ai followed by SU {2) l ® U {1)y , we get the mass terms for the charged leptons and neutrinos. 
Assuming the vacuum alignment 



^) = K,o,o) 



(16) 



the charged lepton mass matrix is given as 




(17) 



Note that we could also obtain a diagonal charged lepton mass matrix even if we assume that e"^, 
and t'^ transform as 1", 1' and 1, and = (0, vt, 0) with appropriate change in the Yukawa 

Lagrangian. Similarly, e^, and transform as 1', 1 and 1", and = (0, 0, Vt) could give us 

the same M.i. In what follows, we will assume that the M.i is of the form given in Eq. (17). 
In the most general case, where all three singlet Higgs as well as 4>s are present and we do not 

assume any particular vacuum alignment, the neutrino mass matrix looks like 



'a + 26i/3 c-63/3 
My = mo\ c - 63/3 d + 262/3 a - bi/3 
d-62/3 a -61/3 + 263/3, 



(1^ 



where mo — ^hi — 2xb^, a — 2xaj^, c — 2x'^^ and d — 2x'^^ and we have written the VEVs as 

(05) = {VS„VS„VS,), iO = U, = U', (D = U", {hu,d) = Vu,d . (19) 



3 Number of Singlet Higgs and their VEVs 

In this section we work under the assumption that the triplet Higgs (f)s has VEVs along the 
direction 



?) = {vs,vs,vs) ■ 



(20) 



This produces the neutrino mass matrix 

M„ = mo 



/a + 26/3 c-6/3 d-b/3\ 

c-6/3 d + 26/3 a -6/3 
^d-b/3 a -6/3 c + 26/3/ 



(21) 



where 6 = 2xb^. In the following we discuss the phenomenology of the different forms of M.u 
possible as we change the number of singlet Higgs or put their VEVs to zero. We assume that 
A^iy is real. 
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3.1 No ^4 Singlet Higgs 

If there were no singlet Higgs, or if the VEV of all three of them were zero, one would get the 
neutrino mass matrix 

/ 26/3 -6/3 -6/3 \ 
Mu^mol -6/3 26/3 -6/3 . (22) 
V-6/3 -6/3 26/3 / 

On diagonalizing this one obtains the eigenvalues 

mi = mo 6, 7712 = 0, ms = mo 6 (23) 



and the mixing matrix 



/ 4 4 ^ 

1 IT — IT 

6 V 3 V 2 



(24) 



Therefore, we can sec that the tribimaximal pattern of the mixing matrix is coming directly from 
the term containing the triplet Higgs (j)s and does not depend on the terms containing the singlet 
Higgs scalars^. However, in the absence of the singlet Higgs contributions to A^,^, the ordering of 
the neutrino masses turn out to be very wrong. In this case, Am|i = —6^ tUq and Am|i = 0, in 
stark disagreement with the oscillation data. 

3.2 Only one ^4 Singlet Higgs 

If we take only one singlet Higgs at a time then there are three possibilities. The resulting 
mass matrices are shown in column 2 of Table 3. One could get exactly the same situation 
with three singlet Higgs and demanding that the VEV of two of them arc zero while that of the 
third is nonzero. The M.^ given in Table 3 can be exactly diagonalized and the eigenvalues and 
eigenvectors arc shown in column 3 and 4 of the Table, respectively. One can see that only the 
case where ^ is present gives rise to a viable mixing matrix, which is exactly tribimaximal [19]. 
One can check that only the first case has the form for M.v given in Eq. (2). Note that each of 
the M.y given in Table 3 possesses an jS'2 symmetry, which reflects the symmetry of Eqs. (10-12)^. 
While the case with ^ exhibits the ^ — r exchange symmetry, the one with ^" remains invariant 
under e — /j, permutation and the one with ^' under e — r permutation. This would necessarily 
demand that while for the first case 6*23 would be maximal and ^13 = 0, for the second and third 
cases ^23 would be either 90° or respectively, and 613 maximal. 

Since only the case with ^ reproduces the correct form for the mixing matrix, we do not discuss 
the remaining two cases any further. This was the model presented by Altarelli and Feruglio in 
their seminal paper [19]. The mass squared differences in this case are 

Amli = (-6^ - 2ab)ml, Amli ^ -4abml , (25) 



^Note that the mixing pattern does not depend explicitly even on the VEV of (/(g. 

^Since we have assumed the vacuum alignment (0s) = {vs,vs,vs), the b terms of Mi, are such that it is 
symmetric under all the three exchange symmetries. 
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Figure 1: Contour plot of Am^i, Am|^ and sum of absolute neutrino masses rrit in the a — b plane, 
for three different values of mo (mo=0.016, mo=0.024 and mo=0.032) for the case with ^ and (j)s- 
The first row shows contour plots for different values of Am^i (in 10~^ eV^), with red dashed 
lines for 6.5 blue solid lines for 7.1, green dashed lines for 7.7, green solid lines for 8.3, and orange 
dotted lines for 8.9. The second row shows contour plots for different values of Amli (in 10~^ 
eV^), with red dashed lines for 1.6, blue solid lines for 2.0, green dashed lines for 2.4, green solid 
lines for 2.8, and orange dotted lines for 3.2. The third row shows contour plots for rrit (in eV), 
with red dashed lines for 0.05, blue solid lines for 0.07, green dashed lines for 0.09, green solid 
lines for 1.1, orange dotted lines for 1.3, turquoise solid lines for 1.4, magenta solid lines for 1.6, 
and brown solid lines for 1.8. 
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Table 3: The mass matrix taking one singlet at a time, its mass eigenvalues, and its mixing matrix. 
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Figure 2: Scatter plot showing regions in a — 6 parameter space for the model considered in [19], 
which are compatible with the current 3cr allowed range of values of the mass squared differences. 
The parameter mo is allowed to vary freely. 
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where mo = ^u/A. Since it is now known at more than 6a C.L. that Aml^ > [7], we have the 
condition that —2ab > b^. Since 6^ is a positive definite quantity the above relation imphes that 
—2ab > , which can happen if and only if sgn{a) ^ sgn{b). Inserting this condition into the 
expression for Am|^ gives us Am|]^ > necessarily in this model. Therefore, inverted neutrino 
mass hierarchy is impossible to get in the framework of the simplest version of the model proposed 
by Altarelli and Feruglio [19]. 

The sum of the absolute neutrino masses, effective mass in neutrinoless double beta decay and 
prediction for tritium beta decay are given respectively as 



In Fig. 1 we show the contours for the observables Am2i, Am|^ and mt in the a — h plane, for 
three different fixed values of tuq. The details of the figure and description of the different lines 
can be found in the caption of the figure. 

In Fig. 2 we present a scatter plot showing the points in the a — b parameter space which are 
compatible with the current 3o" allowed range of the mass squared differences. We have allowed 
tuq to vary freely and taken a projection of all allowed points in the a — b plane. Note that 
while a is related to the VEV of the singlet b is given in terms of the VEVs of the triplet (ps- 
We reiterate the point mentioned before that the TBM form for the mixing matrix comes solely 
from the vacuum ahgnment of (f)s and ^ is not needed for that. The singlet ^ is necessary only for 
producing the correct values of Am'^i and Aml^. However, we can note from Fig. 2 that for a given 
value of a needed to obtain the right mass squared differences, the value of b is almost fixed. In 
fact, we can calculate the relation between a and b by looking at the ratio = — 0.03, 

where 0.03 on the right-hand-side (RHS) is the current experimental value. This gives us the 
relation b ~ —1.88a. Therefore, not only does one need the ahgnment {(ps) = (vs^vs^vs) to get 
TBM, one also needs a particular relation between the product of Yukawa couplings and VEVs of 
(f)S and ^ in order to reproduce the correct phenomenology. This appears to be rather contrived. 
Even if one includes the 3(7 uncertainties on Amli and Am|i, \b\ is fine tuned to |a| within a factor 
of about 10"^. 



If we take two singlet Higgs belonging to two different singlet representations and allow for nonzero 
VEVs for them, then the Aii, obtained for the three possible cases are shown in the first three 
rows of Table 4. One can again see that of the three possible combinations, only the ^" 
combination gives a viable TBM matrix. The other two mass matrices exhibit e — r (^, ^")and 
e — fi symmetry respectively and arc ruled out. Note that we have chosen a = c for the ^, ^" 

combination, a = d for the ^, ^' combination and c = d for the C," combination for the results 
given in Table 4. This is a reasonable assumption to make since the phenomenology of the three 
cases does not change drastically unless the VEVs of the singlet Higgs vary by a huge amount. In 
particular, by changing the relative magnitude of the VEVs, we do not expect the structure of the 
mixing matrix for the first two rows of Table 4 to change so much so that they could be allowed 



mt = \mi\ + \m2\ + Ims], (mee) = mo(a + 26/3), 




(26) 



3.3 Two A4 Singlet Higgs 
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Table 4: Here we take two singlets at a time, and analytically display eigenvalues and eigenvectors 
of the neutrino mass matrix. 
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Figure 3: The left panels show sin^^i2, sin^6'i3 and sin^ 6'23 vs e and the right panels show Am^i, 
Amg]^ and rrit vs e respectively. Here ^' and ^" acquire VEVs. The other parameters c, b and niQ 
are allowed to vary freely. 
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Figure 4: Scatter plot showing the 3a allowed regions for the b—c—d parameters for the case where 
^' and ^" acquire VEVs. The top, middle and lower panels show the allowed points projected on 
the c — b, c — d and d — b plane, respectively. The parameter mo was allowed to take any value. 
Here we have assumed normal hierarchy. 
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Figure 5: Same as in Fig. 4 but for inverted hierarchy. 
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by the current data. In the hmit that c = it is not hard to appreciate that the resultant matrix 
with ^' and C," would exhibit /i — t symmetry, though the ^' and ^" terms alone have e — r and 
e — II symmetry respectively. 

Since the ^" combination is the only one which gives exact TBM mixing in the approximation 
that c = d, we perform a detailed analysis only for this case. Putting c = c? is again contrived and 
would also lead to a certain fixed relation between them and b, as in the only case. This would 
mean additional fine tuning of the paramaters, unless explained by symmetry arguments. Hence, 
we allow the two VEVs to differ from each other so that c — d+e. If e is small we can solve the 
eigenvalue problem keeping only the first order terms in e. The results for this case are shown in 
the final row of the Table 4. The deviation of the mixing angles from their TBM values can be 
seen to be 

C. = 0. D^^-± (27) 

where D12 = sin^^i2 — 1/3 and D23 = sin^^23 ~ 1/2. We show in the left hand panels of Fig. 
3 the mixing angles sin^ ^^12 (upper panel), sin^ ^13 (middle panel) and sin^ 6'23 (lower panel) as a 
function of e. We vary e from large negative to large positive values and solve the exact eigenvalue 
problem numerically allowing the other parameters, mo, b and d, to vary freely. For e = of course 
we get TBM mixing as expected. For very small values of e, the deviation of the mixing angles 
from their TBM values is reproduced well by the approximate expressions given in Eq. (27). For 
large e of course the approximate expressions fail and we have significant deviation from TBM. 
The values of sin^ 9i2 and sin^ 6*13 increase very fast with e, while sin^ 623 decreases with it. We 
have only showed the scatter plots up to the current 3a allowed ranges for the mixing angles. The 
mass dependent observables can be calculated upto first order in e as 

Amli ml{b + d + ^){3d - b + —), Amj^ ~ {4bd + 2be)ml , (28) 

26 v-i , 2 2/26^ ^,0 ^, 26e, 

(mee)=moy, TTit = 2_^\mi\, ~ mo(— + 2d — + 2de-—). (29) 

Of course, epsilon can be larger and we show numerical results for those cases. 



3.3.1 Normal Hierarchy 

The right panels of Fig. 3 show Am^i (upper panel), Am\i (middle panel) and mt (lower panel) 
as a function of e assuming mi < m2 <S rriz- We show Am|i and Am\i only within their Scr 
allowed range. We notice that while Am^i and Am|j are hardly constrained by e, there appears 
to some mild dependence of on it. 

Fig. 4 gives the scatter plots showing the allowed parameter regions for this case. The upper, 
middle and lower panels show the allowed points projected on the 6 — c, d — c and b — d plane, 
respectively. 
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Figure 6: Scatter plot showing the 3a allowed regions in the model parameter space for the case 
where ^, ^' and ^" all acquire VEVs. The upper panels show allowed regions projected onto the 
a — c, a — d, a — b planes. The lower panels show allowed regions projected onto the c — d, c — b, 
d — b planes. Normal hierarchy is assumed. 

3.3.2 Inverted Hierarchy 

For this case it is possible to obtain even inverted hierarchy. We show in Fig. 5 the scatter plots 
showing the allowed parameter regions for inverted hierarchy. We have allowed ttiq to vary freely 
and show the allowed points projected on the c — b, c — d and d — b planes. One can check that 
only for the points appearing in this plot, < mi < 1712. We reiterate that these points also 
satisfy the 3a allowed oscillation parameter ranges given in Table 1. 

3.4 Three ^4 Singlet Higgs 

Finally, we let all three singlet Higgs VEVs contribute to M.^. In this case one has to diagonalize 
the most general mass matrix given in Eq. (21). This matrix has four independent parameters. 



16 



Higgs 



Neutrino mass matrix 



Eigenvalues 



i_ 

V vl 



Eigenvectors 

^73 \ 



1 

1 

V3 



/a + f c-l 
C', mo c-| d+f 



a — 
c+'- 



a — c — d; 

^ mob, \ 

3moa, 
V mob J 



a ^ c = d; 

mo(a + b — c), 
mo{a + 2c), 
m.o{b + c — a) 



a — c^d] 

/ moib + d — a), 

mo (2a + (i), 
\ mo(a + 6 — d) 



a — d^c] 

I moib + c — a) 
mo (2a + c), 
\ mo (a + b — c) 



c — d-\- e] 

I mo{a + b- d - 

mo(a + 2d + e), 
\ mo{b + d - a + I) 



1 




1 

f 









1 






72 




1 


1 


J 


V3 


V2 



/ L A. L \ 



2 J_ 

3 ^3 

1 1 

t/6 a/3 



V2 



75 / 



/ 


i 
1 


i 

1 


i 

1 


\ 






73 






V 


71 


1 
V3 





/ 



+ a2 



-75 + «3 



1 

t 
t 

V3 



a4 



1 

\V2 



as 



72 ~ "6 



Table 5: The mass matrix taking all three singlets, its mass eigenvalues and its mixing matrix. 
The correction factors in the last row are: a2 = , » , as = —-r-^r^^ and Ga = tttst — rr, 

4V2(a-(i) ' 4V2(a-d) * 2V2(a-<i) ' 

"5 = 4V2(a-d) ' "6 = 4V2{a-d) " 
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Figure 7: Same as Fig. 6 but for inverted hierarchy. 

If we assume that the singlet VEVs are such that a = c = c?, ^ then the eigenvalues and mixing 
matrix are given in the first row of Table 5. This gives us a mass matrix whose structure is identical 
to that given in Eq. (22). Hence, it is not surprising that the corresponding mixing matrix that 
we obtain has exact TBM mixing and two of the mass eigenstates are degenerate. Therefore, to 
get the correct mass splitting it is essential that (i) we should have contribution from the singlet 
VEVs and (ii) the contribution from the the three singlets should be different. If we assume that 
a = c 7^ d, then one can easily check that Ai^, has e — r exchange symmetry, and hence the 
resulting mass matrix is disallowed. This is because for a = c, as discussed before we get e — r 
exchange symmetry and the ^' term has an in-built e — r symmetry. Similarly for a = d ^ c, one 
gets e — fi symmetry in A^j, and is hence disfavored. Only when we impose the condition c = d, 
we have fi — t symmetry in Ai^, since the ^ term and the sum of the ^' and ^" terms are now 
separately — r symmetric. Therefore, the case a ^ c = d gives us the TBM matrix and the mass 
eigenvalues are shown in Table 5. 

Since a ^ c = dis the only allowed case for the three singlet Higgs case, we find the eigenvalues 

''Note that one would need additional symmetries to explain a = c = d. 
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and the mixing matrix for the case where c and d are not equal, but differ by e. We take c = d + e 
and for small values of e give the results in the last row of Table 5, keeping just the first order 
terms in e. The deviation from TBM is given as follows 

The mass squared differences are 

Aml^^ml{2a + b + d+^){3d-b+ Aml^ 2mlb{2d - 2a + e) . (31) 

Prom the expression of the mass eigenvalues given in the Table, one can calculate the observables 
rrit, nip and {ruee) 

, . / 2&. , , 2 2 / 2 4a6 26^ ,2 Ahd , 2be. . , 
(mee) =mo(a + y), mt = 2^|mi|, ~ mo (a + — + — + 2c( — + 2de - —) . [32) 

i 

3.4.1 Normal Hierarchy 



Let us begin by restricting the neutrino masses to obey the condition mi < m2 <S m^ and allow a, 
b, c and d to take any random value and find the regions in the a, b, c and d space that give Am2i, 
Amg]^ and the mixing angles within their current 3a allowed ranges. This is done by numerically 
diagonalizing M.^. The results are shown as scatter plots in Fig. 6. To help see the allowed zones 
better, we have projected the allowed points on the a — c, a — d and a — b plane shown in the 
upper panels, and c — d, c — b and d — b plane in the lower panels. There are several things one 
can note about the VEVs and hence the structure of the resultant Mi, 

• a = is allowed, since this gives a, Mu which has contributions from ^' and discussed in 
section 3.3, 

• 6 = is never allowed since b is needed for TBM mixing as pointed out before, 

• a = b, a = c and a = d are never allowed, 

• c = c? is allowed and we can see from the lower left-hand panel how much deviation of c from 
d can be tolerated, 

• c — and d — can also be tolerated when a 7^ 0. 
All these features are consistent with the results of Table 5. 



3.4.2 Inverted Hiercirchy 

In this case too its possible to get inverted hierarchy. The corresponding values of the parameters 
of A^,^ which allow this are shown as scatter plots in Fig. 7. Here mo has been allowed to take 
any value, and we show the points projected on the a — c, a — b, a — d plane in the upper panels and 
c — b, c — d, b — d plane in the lower panels. Each of these points also satisfy the 3a experimental 
bounds of Table 1. Note that for a = we get the same regions in 6, c and d, as in Fig. 5. 
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Figure 8: Variation of sin^ 9i2 with r] for the case where we allow for a misalignment of the triplet 
Higgs such that {(f)s) = {vsi, vs, vs) and a ^ c = d. 

4 Vacuum Alignment of the Triplet Higgs 

In case we do not confine ourselves to {(ps) = iysi^s^vs), we would have the general given 
in Eq. (18). Since we have argued in the previous section that the only viable scenario where 
one allows for all three Higgs singlet is when a 7^ c ~ rf, we will assume that this condition for 
the singlet terms holds. We further realize that to reproduce a mixing matrix with ^13 ~ and 
62^ ~ 45°, it might be desirable to keep fi — t symmetry in the mass matrix. Therefore, we show 
our results for the case {(ps) = {vsi,vs,vs)- The mass matrix is then given as 



Of course for bi = b one would recover the case considered in the previous section and TBM mixing 
would result. The possibility of 61 7^ 6 gives rise to deviation from TBM mixing. In order to solve 
this matrix analytically we assume that bi = b + rj and keep only the first order terms in rj. The 
mass eigenvalues obtained are 



/a + 261/3 d-b/3 d-b/3 
My = mQ d-b/3 d + 2b/3 a-bi/3 
V d-b/3 a -61/3 d + 2b/3 




(33) 



nil = TTT'O {a + b — d -\ — ), m2 = rriQ (a + 2(i), = rriQ (—a + 6 + c? H — ) 



(34) 
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and the mixing matrix is 



id 



3(3d-6) 



2rj_ 



i 1 + 



2r? 



3(3d-6) 
2V 



3(3d-6) 

^; 

3(3d-fe) 
3(3<i-6) 



3(3d-6) 

Therefore, the only deviation from TBM comes in 9i2 and we have 

477 



(35) 



D12 ~ 



9{3d - b) 



(36) 



We show in Fig. 8 variation of sin^ 612 with rj. As expected, sin^ 9i2 is seen to deviate further and 
further from its TBM value of 1/3 as we increase the difference between vs and vsi- The other 
two mixing angles are predicted to be exactly at their TBM values due to the presence of — r 
symmetry in Mu- They would also deviate from TBM once we allow for either vs^ 7^ vs^ or c 7^ d, 
and in the most general case, both. 



5 Conclusions 

Current data seems to be pointing towards the existence of tribimaximal mixing. One needs to 
invoke some symmetry argument in order to get tribimaximal mixing. The discreet symmetry 
group has received a lot of attention in recent times as an attractive option for explaining 
the masses and mixing pattern of the neutrinos along with those of the charged leptons. We 
made a detailed phenomenological study of the viability of the different mass matrices that can 
be generated by spontaneous A4 symmetry breaking. 

In particular, we considered the model proposed in [19] and studied the phenomenological 
implications for it. The authors of [19] consider only one singlet and one triplet Higgs 
contribution to the neutrino mass matrix. Since the singlet transforms as 1 and since they take the 
vacuum alignment {(ps) = (vsjVsjVs) for the triplet, the neutrino mass matrix by construction 
produces tribimaximal mixing. A lot of attention has been paid on justifying the vacuum alignment 
of the triplet Higgs which is absolutely essential for tribimaximal mixing. In this paper we pointed 
out that in addition to the vacuum alignment of the triplet, one also needs a certain fixed relation 
between the product of the VEV and the Yukawa of the singlet and the triplet. In particular, we 
found that in order to generate the correct ratio of Amgi to Am|]^, one demands that b ~ —1.88a, 
where a = 2xa^ and b = 2xb^. This appears to be extremely contrived and therefore undesirable. 
Even if one includes the 3a uncertainties on Am|^ and Am|^, \b\ is fine tuned to \a\ within a 
factor of about 10~^. Even with this fine adjustment of the product of the VEVs and Yukawas, 
one would be able to generate only normal hierarchy for the neutrino mass spectrum. 

We checked if it was possible to generate a viable mass matrix if the singlet Higgs belonged to 
either the 1' or 1" representation. We found them to be unsuitable due to the in-built wrong S2 
symmetry of the neutrino mass matrix for these cases. We studied the possibility of combining 
two singlet Higgs at a time. We showed that the case where ^' and ^" transforming as 1' and 1" are 
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taken together is the only viable option, since this could lead to an approximate jj, — t symmetry 
for the mass matrix. In the limit that c — d,we get exact — r symmetry and exact tribimaximal 
mixing. We allowed the breaking of this n — t symmetry and showed how the mixing angles 
deviate from their tribimaximal values. We gave approximate analytical predictions for ^m\^, 
Amg]^, rrit and m^, when the difference between the model parameters c and d is small. We 
showed numerically through scatter plots how some of these quantities varied, as the difference 
between c and d was increased. We identified the regions of the model parameter space which 
produce values of the neutrino oscillation parameters within their current 3(7 limits. This model 
is capable of producing even inverted hierarchy. We showed the regions of the b — c — d parameter 
space which allow < mi < m2 and at the same time correctly reproduce the neutrino oscillation 
parameters within their 3a ranges. 

We next allowed for all three singlet Higgs contribution to the neutrino mass matrix. The case 
where a ^ c — d emerged as the only possible case which produced exact tribimaximal mixing. 
Allowing c ^ d allows for deviation from tribimaximal mixing. We studied this deviation as a 
function of the difference between c and d. This model can give both normal and inverted hierarchy. 
The regions of the model parameter space which reproduce neutrino oscillations parameters within 
their current 3a allowed ranges were identified both for the normal as well as for the inverted 
hierarchy. 

Finally, we allowed (^5) to deviate from {vs, Vs, Vs)- Changing the vacuum alignment immedi- 
ately affects the tribimaximal character of the neutrino mixing matrix. We showed the results for 
{4>s) = {vsi,vs,vs) and a ^ c = d only, just for the sake of illustration. Since the resultant mass 
matrix by construction still has the residual /i — t symmetry, 623 and ^13 remain at their tribimax- 
imal values of 45° and respectively. We showed the deviation of sin^ ^12 from 1/3, analytically 
for small values of 77 and numerically for all values of rj, where rj quantifies the difference between 
Vs and vs^- The most general case would of course be when vs^ 7^ 7^ ^s^, where all three 
mixing angles will deviate from their tribimaximal values. Also, the condition c = d is extremely 
"unnatural" and any deviation from that would also break tribimaximal mixing. 
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